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We present a 3-D hydro + cascade model, in which viscosity and a realistic freeze- 
out process for the hadronic phase are taken into account. We compare our results to 
^ ! experimental data and discuss the final state interaction effects on physical observables. 
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m : 1. FREEZEOUT PROCESS AND VISCOSITY IN HYDRODYNAMICS 

O _ 

Hydrodynamic models have been very successful in describing the collective behavior 
i/-} , of matter at RHIC, such as single particle spectra and elliptic flow. In particular the 
strong elliptic flow which, for the first time, reaches the hydrodynamic limit at RHIC, 
provides us with a new understanding of the nature of the quark-gluon plasma (QGP) 
created at RHIC as strongly interacting or correlated QGP [[I]. However there exist 
a number of experimental observations that contradict perfect hydrodynamic models: 
transverse momentum spectra above 2 GeV, elliptic flow at large pseudo-rapdities rj and 
Hanbury Brown - Twiss (HBT) interferometry. These observations suggest that there 
exist limitations to the application of a simple perfect hydrodynamic model to RHIC 
physics and that an improvement on a perfect hydrodynamic model is needed in order 
to obtain a comprehensive and unified description of the data from the point of view of 
hydrodynamics. 

In general, hydrodynamic models require initial conditions, an equations of state (EoS) 
and freezeout conditions as parameters besides the relativistic hydrodynamic equation. 
One of the main advantages of the hydrodynamic model lies in its ability to investigate 
the relation between the EoS and physical observables via comparison to experimental 
data. However, it has been pointed out repeatedly, e.g. by Hirano [|2], that details of 
the treatment of the freeze-out process can have large effects on physical observables. For 
example, the data on spectra and elliptic flow show that the assumption of chemical 
equilibrium or partial chemical equilibrium [0] in the freezeout process is not realistic 
[12]. In addition, studies of collective flow at AGS, SPS and RHIC based on both of 
hydrodynamic models and cascade models suggest that the effects of viscosity are not 
negligible [|2]- Therefore we construct a hybrid 3-D hydro + cascade model to include 
a realistic treatment of the freezeout process and viscosity in the hadronic phase. Such 
hybrid models have been implemented in the past, however with reduced dimensionality 
in the hydrodynamic component and thus with restrictions to observables at mid-rapidity 



X 



*Present address: School of Physics and Astronomy, University of Minnesota, Minneapolis, MN 55455. 
tThis work was supported by DOE grants DE-FG02-05ER41367 and DE-FG02-03ER41239. 



2 



C. Nonaka and S. A. Bass 



[ IH] • As a cascade model we use UrQMD in which final state interactions are included 
correctly [|7]. 

2. 3-D HYDRO + CASCADE MODEL 

2.1. 3-D hydrodynamic model 

We solve the relativistic hydrodynamic equation, 

d,T^ = 0, (1) 

where T^ u is energy momentum tensor, including baryon number conservation, 
dn(riB{T, = 0. We modify our original code for the hydrodynamic expansion in 

Cartesian coordinates [El to that in light-cone coordinates (T,x,y,n) (r = y/t 2 — z 2 ), 
in order to optimize our hydrodynamic model for ultra-relativistic high energy heavy 
collisions [ 9J. We employ Lagrangian hydrodynamics which has the following advantages 
over Eulerian hydrodynamics [ 8 a : 1) We can trace the adiabatic path of each volume 
element of fluid on phase diagram. 2) We can easily investigate effects of phase transition 
on physical observables. 

We assume that hydrodynamic expansion starts at r = 0.6 fm. Initial energy density 
and baryon number density are parameterized by 

e(x,V,v) = e mSLX W(x,y;b)H(r]), 
n B (x,y,r]) = n Bm£kX W(x,y;b)H(r]), (2) 

where b and e max (n Bmax ) are the impact parameter and the maximum value of energy 
density (baryon number density), respectively. W(x,y;b) is given by a combination of 
wounded nuclear model and binary collision model [ and H(n) is given by H(n) = 
exp — — ?7o) 2 /2cr^ • 9(\r]\ — r] ) . Parameters e max , n Bmax , n and <t v are listed in Table 
[0 We set initial flow in the longitudinal direction to Vl — T) (Bjorken's solution) and 
vt = in the transverse plane. We use an equation of state with a 1st order phase 
transition, namely a Bag model EoS with and excluded volume correction [ lllj. The 
thermal freezeout temperature is 110 MeV. 

2.2. 3-D hydro + cascade model 

The key improvements of the hybrid model approach compared to the 3-D hydro- 
dynamic model are the inclusion of viscosity in the hadronic phase as well a realistic 
treatment of final state interactions and a species dependent freezeout process. For the 
transition from hydrodynamic model to UrQMD we introduce the switching temperature 
Tsw which should be chosen below but near the critical temperature (T c = 160 MeV [Ej). 
We set the switching temperature to 150 MeV. The calculated procedure is as follows: 
first hadron distributions from the hydrodynamic model are calculated at the switching 
temperature via the Cooper- Frye formula [Q21- Second, initial conditions for UrQMD on 
an event by event basis are generated from the hadron distributions via the hydrodynamic 
expansion through a Monte Carlo implementation. Finally the UrQMD calculation starts 
with these initial conditions. In this scheme we neglect the reverse contribution from 
UrQMD to the hydrodynamic model. For an initial condition for 3-D hydro + UrQMD 
we assume the same parameterization as that of pure hydrodynamic model except the 
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value of maximum energy density, which we have adjusted to achieve a better agreement 
to experimental data. The parameters are summarized in Table 1. 

Table 1 

Parameters for initial conditions of pure hydro and hybrid model. 





r (fm/c) 


max. of e (GeV/fm 3 ) 


max. of fi Bmax (fnr 3 ) 


Vo 


a v 


pure hydro 


0.6 


55 


0.15 


0.5 


1.5 


hydro + UrQMD 


0.6 


43 


0.15 


0.5 


1.5 



3. RESULTS and DISCUSSIONS 

Figure [T] shows the Pt spectra of ir, K, and p at Au + Au y/s = 200 GeV central 
collisions based on 3-D hydro + UrQMD. Our calculated results are very close to PHENIX 
data [112] up to P T ~ 2 GeV. In 3-D hydro + UrQMD model the hadron ratio 7r/p and n/K 
are obtained correctly. In the case of the pure 3-D hydrodynamic model with low thermal 
freezeout temperature under the assumption of chemical equilibrium, we fail to obtain 
the correct hadron ratios This observation emphasizes the importance of the realistic 
treatment for the freezeout process, but further studies are necessary to investigate the 
effects of cross-over transition as well. 

Figure El shows the centrality dependence of Pt spectra of 7r + . We can see that in 
peripheral collisions the difference between experimental results and calculated results 
appears at lower Pt compared to central collisions. These deviations are indicative of 
the diminished importance of the soft, collective, physics described by hydrodynamics 
compared to the contribution of jet-physics in these peripheral events. 



."10° 
^ 10 2 

o 

£ 1 
I 10 ' 1 

5 10* 

gio" 3 

^10" 4 



• 


centrality 0-6 % 
O 7r + (PHENIX) 
A K + (PHENIX) 
□ p (PHENIX) 
• ir + (Hydro + UrQMD) 
A K + (Hydro + UrQMD) 
|. ■ p (Hydro + UrQMD) 




% . 



- 10° 

i 10 2 

£ 10 

■o 1 



% 10 



iio- 2 
£io- 3 
4io- 4 







1 



10" 



••i •••• 



O PHENIX 0-5% 
O PHENIX 10-15 % 
PHENIX 15-20 % 
O PHENIX 20-30% 

• Hydro + UrQMD 

• Hydro + UrQMD 
Hydro + UrQMD 

• Hydro + UrQMD 

' S8 iS 

•s. « as 



See. 1 !?! 







1 



2 

Pt [GeV] 



Figure 2. Centrality dependence of Pt 
spectra of 7r + with PHENIX data 



Pt [GeV] 

Figure 1. Pt spectra for tt + , K + , p at 
central collisions with PHENIX data [ I13j . 

Figure El shows the centrality dependence of the pseudorapidity distribution of charged 
hadrons together with PHOBOS data [E]- The impact parameters are set to b = 
2.4,4.5,6.3 and 7.9 fm for 0-6 %, 6-15 %, 15-25 % and 25-35 %, respectively. Our results 
agree with the experimental data not only at mid rapidity but also at large rapidity. 

Figure 13] shows the mean Pt as a function of rapidity for n, K and p. The difference 
between open symbols (hadronic decays without rescattering) and solid symbols (full 
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Figure 3. Centrality dependence of rapid- 
ity distribution with PHOBOS data 



0.5 



0.0 



• i burner A BRAHMS 0-5 % n 
O T + decay A BRAHMS 0-5 % K + 

• K burner 
O K + decay 

p burner 

O pdecay » m ****** 



-4 -3 -2 -1 





y 



12 3 4 



Figure 4. Mean for tt + , K + , p as a 
function of y with BRAHMS data [IIS]. 

hadronic rescattering) determines the reaction dynamics of the final state interactions. 
We find that the average transverse momentum (Pt) for 7r decreases whereas the proton 
Pt increases. Since the number of ir is larger compared to that of p, the effect on the 
protons per particle is larger than for the pions. This transfer of transverse momentum 
is commonly referred to as the " pion-wind" . 

In summary, we have presented a fully 3-D implementation of a hybrid hydro+micro 
transport model. The model includes a realistic treatment of hadronic freeze-out as well 
as an implementation of viscous effects in the hadronic phase. We have compared single- 
particle spectra and pseudo-rapidity distributions to data and have studied the reaction 
dynamics of protons, pions and kaons in the hadronic phase of heavy-ion collisions at 
RHIC. 
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